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ABSTRACT 


The instantaneous force acting on a circular cylinder 
immersed in a unidirectional periodic flow characterized by 
U=VFt wh Sin 2mft was investigated. The total force was 
divided into three components consisting of a steady-state 
drag; an oscillating drag which is a function of the 
frequency parameter A = fd/V ; and an inertial force which 
is also a function of the frequency parameter. The three 
resistance coefficients associated with the force components, 
Cy(0), CHA), and Cy 6A) were calculated. 

It Was Peulideenate the resistance coefficients» aremnet 
materially affected by the free-stream oscillations whenever 
the frequency of the oscillation is lower than the natural 
frequency of the vortex shedding. However, when A approaches 
the Strounal Prequencyeot vortexeshedd ineyethe= inertia 
coetiricilene Cu is significantly reduced, and likewise when 
X approaches the value of twice the Strouhal number, CpOA) 
is markedly reduced. inereswoms below unis “range, “the 
inertial force becomes negligible and the total resistance 
could be Aiewented accurately by the coefficients Cy (0) 


and CHOA) with C,, fixed at its theoretical unseparated-flow 
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Symbol 


C0) 
Cy) 


NOMENCLATURE 


effective cylinder flow-blockage area (A = ak) 
instantaneous acceleration 

cylinder cross-sectional area (A, = nd</4) 
cross-sectional area of test section (A = ht) 


ES 


steady state drag coefficient 

coefficient for steady state portion of drag 
CocmmencCn lL enor Ose) lating portion of drag 
COcmmvement for inertia force 

Giameter of test cylinder 

wire Lecorrecuson factor 


instantaneous total force acting on the test 
cylinder 


frequency of oscillation 
dGmac forces: acting on the test cylinder 
inertia forces acting on the test cylinder 


hemeht of-tunnel test section 


constant (K, = %V°A) 
constant (K,, = VV Ap) 
2 
= lk 
constant (K, aVy Ae) 
constant (K), = emfpA &V,) 
width of tunnel test section; test cylinder length 


frequency of vortex shedding 
pressure 


height of dynamic pressure curve at t = T/4 





Re 


Re 


St 


<= <= 4 c¢ 


K Mm <j 


measured velocity 
test cylinder radius 


height of dynamic pressure curve at t = 3T/4 


Reynolds number (Re = Ud/yv) 
Reynolds number (Re = Vd/v) 
Strouhal number (St = nd/V) 


permed ofeesed | lataen 

time 

instantaneous velocity (U = V + Vy sin27ft) 
free stream velocity 

amplitude of velocity oscillation 

average velocity 


horizontal distance between pressure taps 1 and 3 


horizontal distance between pressure tap 3 and 
center of test cylinder 


normalized frequency of oscillation (A = fd/V) 
fitted Kinematic viscesmty 


fluid density 





Ii. IMGRODUCTION 


A. EXPERIMENTAL JUSTIFICATION 

The subject of forces acting on bluff bodies immersed 
in time-dependent flows has been aa will continue to be of 
gnverest to fluid dynamicists, aerodynam#eaists, and practicing 
engineers for the special reason that in mature neither the 
body nor the fluid which surrounds the bo@g is ever in a 
state of steady motion. Even the flow behiind a bluff body 
moving steadily through a fluid is accompanied by large 
scale unsteadiness. Thus, any type of unsteadiness of the 
ambient flow and/or the motion of the body introduces 
eae tional changes in the characteristics of the Ilow and 
its analysis. 

The drag and inertial forces are jemeependent as well 
ao time dependent, and the resistance coed@icients obtained 
in unseparated flows are not applicable to separated cases. 
Although indirect, the role of viscosity £8 paramount in 
that its conte Gites are separation, vortex formation and 
SAeacing, and resultant alternations in the virtual mass. 
It is thus clear that it is necessary to determine the 
relationships between various resistance components in terms 
of the unsteadiness of the ambient flow and/or body motion, 
Beometry of the body, the degree of upstream turbulence, 


meet history of the flow, etc. [Ref. 1]. 





However, the understanding of the behavior of time- 
G@ependent fiows about bluff bodies may be advanced only by 
Sensiaering the relatively more manageable cases and gradually 
integrating the information so obtained. This led to the 
present experimental investigation of the unidirectional 


meriodic flow about circular cylinders. 


Ee SURVEY OF THE PREVIOUS INVESTIGATIONS* 


A review of some of the previous investigations on tTime- 
G@ependent flows) may best be presented by dividing the motions 
under consideration into various but admittedly arbitrary 
categories. 


i. Theoretical Analysis of Unseparated Time-Dependent 
Laminar Flows 


The study of the theory of separation-free time- 
dependent laminar flows has enjoyed particular attention, 
Pertiy dué to its practical significance and partly due To 
its relative mathematical simplicity. Surveys by Stewartson 
{Ref. 2], Stuart [Ref. 3], and Rott [Ref. 4] reflect the 
current level of understanding. 

The unseparated class of unsteady flows, most of 
Pench results fren the unidirectional or periodic acceleration 
of bodies in an infinite or bounded fluid medium, gives rise 
to an induced mass which must be added to the real mass of 


ome body. Inweeneral, the virtual inertia tensor depends 


This section has been reproduced here, by permission, 
from an unpublished report by Prof. T. Sarpkaya. 





@mene Shape of the body, the nature of the fluid mediun, 
wae orientation of the moving ObmeeL With respect to the 
direction of motion, and the depth of relative submergence 
from a free surface and/or solid boundary. 

The effects of viscosity and separation on added 
mass or the dependence of added mass of a given body on the 
characteristics of the time-dependent separated flow cannot 
wee be theoretically evaluated nor can it be experimentally 


separated from the total resistance experienced by the body. 


oO. ite lL ilaimplitude Vabration of Bodies in a Liquid 


Otherwise at Rest and the Initial Instants of 
Acceleration of Bodies from a State of Uniform 


Velocity 


Vibratory motion has often been used [Refs. 5-8] as 
one of the experimental means to determine the virtual mass 
coefficients. The technique has been restricted to amplitudes 
of motion which are so small that either separation does not 
occur or the characteristics of separation do not change 
due to the imposed, high-frequency vibration. Thus, the 
results obtained are not applicable to occurrences in which 
eae Guravion ot unsteady flow in one direction is long enough 
for separation to occur. 

It is not necessary that the body be initially at 
rest to determine the added mass provided that one starts 
with a steady ambient flow about the body or with the body 
moving at constant velocity in a fluid otherwise at rest. 
Hamilton and Lindell [Ref. 9] have shown that the added-mass 


coefficient of a sphere determined by imposing a uniform 


10 





peeelecravHon ohito ies» steadyestate motion #s equal to its 
theoretically determined value as well as to its experimen- 
tally determined value obtained by small-amplitude, high- 
frequency oscillations. 

Rayleigh [Ref. 10] has shown that if the acceleration 
continues, i.e. V = V + kt”, both the drag and the added-mass 
eoeffierzents change with time. 

The foregoing arguments suggest that the force 
exerted by the fluid on a body immersed in a time-dependent 


flow might be expressed as 


F = 4CypA|V|V ct Cy, av/dt + history-dependent (1) 


drag and inertia forces 


in which A is the projected area of the body and M the mass 


of the fluid displaced by the body. The drag coefficient Ch 
is to be taken equal to the quasi-steady state resistance 
coefficient and Cu to be taken equal to that measured or 
calculated from the potential theory when the value of the 


history-dependent force is zero [Ref. 9]. 


3. Bodies Subjected to Large Amplitude Harmonic 


Oscillations in a Fluid Otherwise at Rest 

A number of significant studies [Refs. 11-15] were 
made to determine the forces acting on bluff bodies under- 
going harmonic, large-amplitude oscillations in a liquid 


otherwise at rest. 


10aR 





Odar and Hamilton [Ref. 11] measured the force acting 
on a sphere subjected to harmonic oscillations and proposed 
a force equation comprised of mhmec sperts: a steady drag 
Gempeonent where the drag coefficient is the well-established 
coefficient for steady translation of a sphere; an inertial 
Seemponent whose coefficient is determined experimentally as 
e“Punetien of v* /ad (a: acceleration; d: diameter of sphere); 
ana a history-dependent component whose coefficient also is 
eetunerion of V°/ad. 

4. Bodies Subjected to Wave Motion with Zero Mean Velocity 

Considerable work [Refs. 16-21] has been done on this 
type of time-dependent flow in an attempt to predict the wave 
hOrces: @ia je lees and other submenmed ssieeuctbusees., dn eneral, 
Ghe tetal force acting on the body is assumed to be composed 
of a velocity-dependent and an acceleration-dependent force. 
The resulting equation, known as the Morison equation [Ref. 


A 


16], is given by 


=aabe Val eave 

F aC pA|VIV + Cym, dV/dt 2) 
where Cur = let Cy Cebeimeg the added mass coefficient. Experi- 
ments show that Cy and Cu depend on time and show considerable 


scatter about their assumed mean values [Refs. 17, 20] over 


ercmele. 


lee 





o Someoon Ubscevreo tO Urmdi rectional Acceleration in 


a Fluid Otherwise at Rest or Unidirectional Unsteady 
Flow About Bodies Held at Rest 


Upstate weeett orl wide Deen made FCO determine the 
components of force acting on a body accelerating unidirec- 
tionally in a fluid otherwise at rest. Iverson and Balent 
[Ref. 22] and Keim [Ref. 23] towed spheres, disks, and cylin- 
ders in unsteady motion through still water. It was concluded 
Chat the total instantaneous force on these objects could be 
Gescribed by a single coefficient of the form C = 2F/pAVe 
and that this coefficient was a function only of a so-called 
meceleration modulus, aa/V°. Keim, however, detected an 
Meparent Reynolds number effect for cylinders. 

Laird et. al. [Ref. 24] measured the forces acting 
on accelerating cylinders and found strong evidence of devia- 
mrorn Of Lhe drag coefficient from the accepted values for 
fromm moulon. They further found that the acceleration 
meauilus did not correlate the résistance coefficient near 
boundary-layer transition. 

Sarpkaya and Garrison [Ref. 22] working with uni- 
Serectional.flow with constant acceleration about cylinders 
and plates have found that both Cy and Cu 


history of motion, that the Morison equation (Eqn. 2) is 


depend on the 


Polidsonly for flows with constant acceleration, and that 


the acceleration modulus can correlate the data only if the 


acceleration is kept constant. 


ays 
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6. Forced or Self-—Induced Transverse Oscillations of 
aepoayeinm a Fluid in Steady Motion 


Experiments have been directed to the effects that 
attend the oscillations of a body in the plane of the lift 
force. This type of interaction becomes particularly impor- 
tant when the oscillation of the structure is the necessary 
Comet On Loretie femeracvcion of the excititas forces. The 
most pertinent information from these experiments is that 
the lift and drag forces act on the oscillating cylinder at 
the vortex shedding frequency and twice the shedding fre- 
G@uency, respectively, provided that the driwing frequency on 
the cylinder is appreciably different from the shedding fre- 
quency. When the forcing frequency of the cylinder approaches 
mae womeagdime frequency, the natural sheddimg frequency is 
lost and it "locks-in" to the forcing frequency. This 
synchronization persists over a range of frequencies which 


may be termed the "range of synchronization" [Ref. 26]. 


ie Uadediierect ional Oscillatory Flow About Bodies at Rest 
and Steady Uniform Flow About Bodies Subjected to 


pureamwise Oscillations 

Relatively re studies have been conducted in this 
category. Chen and Ballengee [Ref. 27] examined the vortex 
Piecdaane from circular cylinders in an oscillating free 
stream. Their results suggested that "in an oscillatory 
freestream of 3 HZ and Reynolds number up to 40,000, the 
vorvex shedding from a circular cylinder responds” instantan- 
eously to the freestream variations" and that "the instan- 


meneows Strovlal nufiber stays sensibly constant at 0.2 + 0.01." 


14 





Hatfield and Morkovin [Ref. 28] studied the effect 
of an oscillating free stream on the unsteady pressure on a 
eoreulear Ccylimder. Mey have found that @here is no signi- 
ficant coupling between the small-amplitude free-stream 
eee | latdons and theevertex shedding.» Tredr results-would 
Suggest that the drag coefficient associated with the mean 
flow would essentially remain constant at its steady state 
value. On the other hand, Mercier [Ref. 29] who subjected 
cylinders to large streamwise oscillations found that the 
average drag coefficient significantly increases with fd/V 
and that the rate of increase depends on the amplitude to 
Gzameter ratio. Combined, these results would suggest that 
the deareemed couplime bebween aoe err Of sesciml lation 
ema the vortex—-shedding frequency and hemee the forces acting 
on the body strongly depend on the amplitude of oscillation. 

Davenport [Ref. 30] subjected bluff bodies to small 
amplitude oscillations in a water flume and evaluated the 
Gvas and inertia coefficients through the @easurement of the 
rate of damping of the amplitude of oscillations of the 
edt a6. i.e. without measuring the forces acting on the 
bodies. Using the frequency parameter fd/V to correlate 


the data, he found large variations in Cy and Care 


ae oCOP RO THE PRESENT STUDY 
The present investigation is limited to an experimental 
determination of the instantaneous forces acting on a 


circular cylinder immersed in a unidirectional periodic 


an 





Plow characterized by 


U= V+ V, sin 2nft. (3) 

Its theoretical aspects are structured along the lines of 
the unpublished work by Sarpkaya [Ref. 31] conducted at 
Aerodynamische Versuchsanstalt~-Gottingen with respect to 
Pie SUDDOSICT1OnN that the instantaneous foree acting on an 
Soject immersed in a periodic flow strongly depends on the 
parameter fd/V and to varying degrees on Vd/v and V/V 
where f is the frequency of oscillation and d the cylinder 
diameter. | 

(ie=ranze ©f these parameters during the experimental 


work was: 


O.@@il < fd/V < 0.05 
LIeOOO< VEY y < "66.000 
0.085 < V,/V < 0.69 


As presented in Section IV, the parameter fd/V is extended 


to a value of i Susanne vnewdata of others, 
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II. APPARATUS AND PROCEDURE 


A. APPARATUS 
1. Water Tunnel and Oscillating Mechanism 

The experiments were conducted in a 500-gallon- 
eapacity, recirculating-type water tunnel as shown in Figure 
1. Water was circulated by a low-rpm, high-capacity, centrif- 
Weal pumpwwnich Weeemecapable of providing test-section mean 
velocities of 2-25 feet per second. The test section 
consisted of a rectangular conduit 4 inches wide, 8 inches 
hieh, and 16 inches long. 

Periodic variations in water velocity were provided 
by a rotating butterfly vaive driven by a variable speed 
motor capable of a range of 0-40 — (Fie. 2). The “maga 
tude of the velocity fluctuations was dependent upon the 
mes in errective head through the yalve. AS this is partially 
a function of the valve-to-conduit area ratio [Refs. 35,36], 
an access plate was placed on the tunnel in order to change 
the effective valve area by adding concentric rings to the 
valve perimeter without removing or changing the valve body. 

Associated with the tunnel is a 15-gallon head tank 
located le feet above the test section centerline. 

e, lest Specimens 

The test specimens consisted of polished aluminum 

cylinders placed across the width of the test section perpen- 


ereular to the direetion of flew (Fig. 3). A small bore 


até 
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hole drilled radially from the cylinder surface at the mid- 
feeemoneservea as the total pressure tap. The test cylinder 
was mounted with SA. aii by a self-—-al@ening bearing fit 
mymeone plexiglass window, with the other end passing through 
the other window to a cantilever beam force-measurement 
system. Two cylinders of diameters 0.75 and 1.0 inches were 
weed. 
3. Instrumentation 

The total force acting on the cylinder within the 
test section was measured directly by a strain-gage fitted 
cantilever beam (Fig. 4). The bridge circuit consisted of 
peur active strain gages, with the 6utput signal being fed 
me a dwel-cheanme!l strip chart recorder. 

Three pressure taps were utilized for velocity and 
aeeelemation caiculations. Two variable-sensitivity, 
Hii een ene pressive transducers connected to two dual- 
channel strip chart recorders were used for this purpose. 
Each transducer/recorder system was calibrated with a water- 


manometer. 


B. PROCEDURE 
1. Measurement of erce 
h COmvImUCUSsUracemol syncs voual ,Sanstantancous force 
acting on the cylinder was obtained from the strain-gage 
fitted cantilever beam (Fig. 4) which was wired to a strip 
chart recorder. This system was calibrated by static loading 


of the cylinder. Three different cantilever beams were used 


Zl. 








SIDE VIEW 


Figure 4: 
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during the course of the investigation, each beam being 
several times more sensitive than the previous one, in order 
to increase the range of experimentation. 
2. Measurement of Frequencies 
Helen Cuerpo welars recorder Was tne Capability of 
variable, controlled paper speeds. Thus, with aac speed 
and distance between corresponding points on the trace 
Pelevni, vie period of the cyclic signal could be determined. 
The frequency of oscillation was then simply the reciprocal 
of the period. 
3. Measurement of Acceleration 
Acceleration, with the exception of a scale factor, 
was measured directly from the differential pressure taken 
from two static pressure taps placed a known distance apart. 


Peetecuing frieccvion losses, the equation of motion yields 


Orie. 5), 
Che P2) Ang = ma (4) 


where m is the mass of fluid. Writing m= po X Ams aig CE 


one has 


a 5) 
or 


Voe (6) 


a5 





4, Measurement of Velocity 


The dynamic pressure, in the form Ween was 
Gevermined using the differential pressure from the total 
and static pressure taps 2 and 3, respectively, (Fig. 5), 
Pyeproperily locating the wall pressure tap for a given size 
ey linder . 

Determination of this distance begins with the 
unsteady Bernoulli equation 

DV 0V 0V 1 oP (7) 


De 9t “ds ~ ~ 9 3s 


which reduces to 


= 0 V 
Upon integration, (8) becomes 
} 
Ape IV 
(CP tape / 2) = -p J a7 ds (9) 
Ab ale 
or 
Poet ov = fp ei ema). 2/2 = - y av ds (10) 
oa te Tc aaa es Obes 


ies eCvacenu from Figure 5 that V5 = 0 at the stagnation 


point. Therefore, (10) becomes 


@ 


S.. 2 2 2 wee eamee oo 
1 


ae det 0M 2 Gis 
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Now, from Equation (6) 
_ an oV 
AP = Pp Ee = (0s aE (12) 


emo combaning these pressure differences yields, 


2 
rr: av av 
P, - Ps = oV /2 - _ ae as + pX ye: Ges) 
In order to satisfy the requirement that P, ~ P. ~ pv" /2. 
it is necessary, considering (13), that 
a a (14) 
mG ot . 


Hor ecvelwuate thas inGegral, the variation of V with s along 
the stagnation streamline must be known. 

The complex potential for a cylinder in a uniform 
stream is given by : 


ae. 
C= WCU + 5) G5) 


and the complex velocity is 


2 
_ A _ a 
i =o WL = 3) | (16) 
But, by definition 
anit Cie 


26 





Tmeretore, "along the Stagration streamline connecting the 


points l and 2, one has 


v= 0 
Ro 
= V(l1 - =) (18) 
Ss 
eile = —Seee USUI ae Enis result into the integral 


ibaa Oeibes a lids 


2 2 2 
if; av ds = CA if (ai = R) ds. Cale 
ot 3 2 
Ik ab Ss 


Camrying out the integration and equating the results with 


the requirement of (14) yields 


2 
2 eet: 7 em acc 

te Stee 2 t SP 1 wee) 

which simplifies to, 
Ro Ro 

Fe Oe Oar ee (2m) 

OF 
Ro 


where (X + Y) is the distance to the point 2 and R is the 


meadius of the cylinder @nd distance to the point l. 
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Since the flow upstream from the cylinder extends 
to relatively large distances, i.e., as X + ~, Equation (22) 


reduces to 
Y = 2R. (23) 


Therefore, pressure tap 3 (Fig. 5) must be located one 
diameter upstream from the center of the cylinder in order 


e 
we ODtLaIN P., - P. = py /e. 


3 

The length of the test section upstream from the 
cylinder was sufficiently large to ignore the small influence 
oe CoOmerdeupon Of “Chewitowmedr the test Section on the 
amerent TfTow. it is Tor this reason that iss taken to be 
eiIfectively infinitely large. 

>. Hydrodynamic Tunnel Correctiom Factors 

oteady state tunnel correction factors as discussed 
jn Pope (Ref. 34] account for the fact that the bluff body 
is not being tested in an infinite medium. There is a 
blockage-effect correction to the measured velocity, which 


ioe Cuan atieets the calculated drag coefficient. In a usable 


fomm, this correction factor is given by 
(24) 


where d is the cylinder diameter and h is the height of the 


tunnel test section. This correction factor is then used to 
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modify the measured velocity, q, by 
U=q{lte) — (25) 
which then modifies the uncorrected drag coefficient to 


C = 2F/pU“La . (26) 
Gomi 
tee Cifect of ths correction is to increase the measured 
PerrOCiLY@and Umais CO @Gecregse the drag coef#ricient. 

However, with the critical placement of the static 
paesoure Lap Seand its proximity to the test’ cylinder, 
another correction must be applied to the measured velocity. 
From Figure 6, for steady state flows, the ratio of local 


fluid velocity to the free stream velocity is given by 


ra 


U = | ——— 2 -___—_, (27) 
R R 
= ote COS SE 7, 
: r Ae 


For the geometry in this situation, the effect of this 
correction is to decrease the measured velocity and thus 
to increase the drag coefficient. 

As the effects of (25) and (27) counteract each other 
in the steady state flow situation, and since no information 
is available for time-dependent flows, neither correction 
was applied. This is particularly justified in the present 


study since the time-dependent drag coefficients were compared 


aS 
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(eeu Calculated Steady drag coefficients, neither being 


Corrected for sany hy@m@eeynamic constraints. 
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III. THEORETICAL CONSIDERATIONS 


\ 


Assuming that the unidirectional periodic flow in which 
the test cylinder is immersed is a sinusoidal flow superposed 
Pamoeoucdey ilow, the instantaneous free stream velocity can 


be represented as 

Ue ee pein en 66. (28) 
iieeretore, the instamtaneous acceleration becomes 

a= anfV , cos 2n ft (29) 


and the velocity squared is 


Uc = 7 + 2, Sapeeen = fb + Va" eine anft. (30) 


Now, in dealing with the total force acting on the 
cylinder in time-dependent flows, it has been most advan- 
Geaseoule vO COnsiger it as composed of a drag force and an 


i”ervciat rorce, such as 
Nl 7 (31) 


Then using the concept of the normalized drag and inertia 


mesmegames COGmriclents, C. and C respectively, the total 


D Mie 
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jiemeewcan be represented as : 


Zz au 


+ C, L0A == (32) 
iM 


= 1 
F C.4ApU Mie 


D 
where A is the effective cylinder flow-blockage area, ne 
is the cylinder cross-sectional area, 2 is the cylinder 
length, and p the fluid density. Upon substitution of (29) 


and (30) into (32), one obtains, 


sin -onft] 


(33) 


F = CLA V + ApVV.sin2nrft + %ApV 


A A 


+ CylofA 2nfV,cos2nft]. 


A 
bewever, tne results of the works previously cited indicate 
UMeac a two parameter total-force equation does not seem To 
be capable of providing sufficient generality to cover the 
conditions encountered in arbitrary time-dependent flows. 
Also, the past history of the flow must be considered along 
with the velocity, rate of change of velocity, and Reynolds 
number. Therefore, defining C,(0) as the steady state 
coefficient of drag which occurs at ane t = 0 When U = V, 
and assuming that a normalized frequency parameter i = fa/V 
provides the necessary correlation between the resistance 
coefficients and the aforementioned flow parameters, the 


total force can be written as 
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rahe C.(0)%ApV 


+ C(A)LVW,Apsinantt + %ApV entre) (C34) 


A 


~ Cy OAILA plentV eeseait,). 


A 
The coefficients C(0), CHOA); and Cry A) are determined 
ceperimenvally. it is obvious that the validity of the 
feregoing analysis depends upon the existance of a correla- 
mom beeween Une Said coefficients and the normalized 
@sciiietion perameten, A, for various walues of Reynolds 


number, A, and V/V. 
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IV. “SVU ATION @ND PRESENTATION OF DATA 


A. RAW DATA 
As cited earlier, the differential pressure directly 
provides the velocity in squared form. Figure 7 shows a 


theoretical curve of ue = (V+vV Sin2nft)© vase Y = 1.5 


A 


and Vy = 0.5. The distances labeled p and r occur at 


t = T/4 and t = 3T/4, respectively. With these distances 


known, V, V and tie at time t = T/8 can be calculated 


ss 
foam Olvouea.  Wraitamaeg at the specific times, 


t = T/4 uv = (V + al = p C359) 

t= 37/4 Uo = (F - v,)° = 7 (36) 

t= 7/8 vue = (¥ + Vv, v2/2)° (37) 
one has 

a ee Vp (38) 
and 

V-V,= ve (39) 
which yields 

V= (Vp + yYr)/2. (40) 


3p 
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af h/Le c/iL h/b S/iL 0 


ee ee 
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Now, upon subtracting (36) from (35), one obtains 


ee oe v,)° ni = v,)° = 4Wv, (41) 


which reduces to 


V, = <n r)/4V. (42) 
Also upon expansion of equation (37), at time t = T/8, one 
has 


e 


Ue oy 2 


“+ bi, + ¥V, (43) 
Figures 8 and 9 display copies of strip chart recordings 
Gime Di. CcCalmdava Seus tor low and high @requencies of os¢cilia-— 
tion, respectively. The upper traces represent 
AP5_2 = pU“/2 while the lower traces indicate the instantan- 
eous total force... Due to the fact that the torque curve of 
a butterfly valve is not symmetrical about 45° nor the 
discharge coefficient (and thus the head loss) symmetrical 
about 90° in unidirectional flow [Ref. 33], the velocity 
curve is not a pure sinusoid. Therefore, the following 
procedure was used in reducing the raw data to usable form: 
1. The period, T, was measured between any two corres- 
ponding points on the curve. Thus T was the true period of 


UReEwreperioive silpmal . 
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Pee emer stahicesmimana r were measured to the highest 

eme lowest points on the pressure curve, respectively. 
The occurrance of these points were labeled T/4 and 3T/4 
iiomeapective Ol theamr distance»on the time axis since the 
minamum and maximum velocities will oceur at these times. 
be | Che hetehiwas the sduianed—-velociwpy curve at times 
t = 0 and t = 7T/8-were cakeulated from p and r as if the 
VMemecity were atime sinusoid. Tae corresponding points on 
Maem curve were then found at these heights, and the base 
line marked 0 and T/8. 

4. Vertical lines were then drawn from the pressure 
curve to the force curve at the places marked 0, T/8, and 
T/4, and the height of the force curve tabulated at these 
points. Since the two traces are made simultaneously and 
m=Oresent INStaMeancous Valwes, these quantities then give 
the forces which would have occurred at times 0, T/8, and 
T/4 had the velocity been a true sinusoid. 

Daeeliemyciwes Of vellocdiymandiforceamere then calcu- 
Herca 1 fom one Calibration, cquations. Reynolds number 
or on the mean velocity was also calculated. For this 
purpose, an average water temperature of 70° F was assumed 
1Or evaluating the density and viscosity of the tap water 


in the tunnel. 


B. EVALUATION OF COEFFICIENTS 
hs dabei Wed in Part IIL, the total force acting on the 


cylinder is assumed to take the form 
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tx 
it 


Cy (0) [¥ApV* 


mea ee | (44) 


+ 


CyoOICVW,Apsinanft + %ApV, 


+ Cy (A) L2mA p£fV, cosenft | 


After the reduction of the raw data, at times t = 0, t = T/8, 
and t = T/4 all the variables with the exception of Cy (0), 
CHO); and Cy OA) are known. However, a system of three 
equations for the three unknowns can be developed. Re- 


writing (44) as 


rx] 
| 


= Cy (0)K, 


sin°2nft } (45) 


oe 


CyOA)[K,sinenft 1 Ke 


+. 


Cy (AI LK, cosenft J 
emereValuating (415) at times 0, T/4, and 17/8 results in 


t= 0 = Cy (0)K, + Cy CA)Ky (46) 


ct 
i} 


ES Reece CL COOIIK: V2/2 + K%] (47) 


+ Cy (A)Ky V2/2 


tan / Ln = C(0)K, + CyOA) EK, + K (48) 
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These three equations can be solved simultaneously for Cy (0), 
Cy(A), and Cy 6A). 

However, in the range of experimentation of this 
Eeuay, the ravio o6f inertia forces to that of the drag 
memees was << |. this fact led to problems of instability 
in the matrix used to solve for the inertia coefficient. 
For this reason, Cay 6A) was fixed at a chosen value, and the 
values of C(O) and CX) were solved using the instantaneous 


data at times t = 0 and t = T/4. The two resulting equations 


were 
FE’ -sarCoek 
_ 7 M4 
t= 0 o)(0) = —_~MA (49) 
1 
and 
F - C(0)K, 
t = T/4 Cy (A) = (Ka ¥ Ka) (50) 


The coefficients Cp (0) and C6) were calculated from 
(49) and (50) by using two different values for Cy 6A). 
First, a theoretical value of Cu 6A) = 2 for unseparated 
potential flow and second, a semi-theoretical and semi- 
experimental value of Cry (AD = 1.5 (Refs. 25,35] were used. 
The resulting values of C(0) and CA) were essentially 
unchanged. The difference between C(O) with Cy 6A) = 2 
and C,(0) with Cy 6A) = 1.5 was negligible to the second 


decimal. This also held true for CfA). 
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C. NORMALIZED PLOTS 

The non-dimensional drag coefficient C(0) is pre- 
sented in Figure 10 as a function of the Reynolds number 
based on average velocities. 

The coefficients C(O); CyOA); and Cy 6A) are then 
presented as functions of the normalized frequency parameter, 


X (Figs. 11-15). 
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V. DISCUSSION OF RESULTS 


A. UNIDIRECTIONAL PERIODIC MOTION 

Mivmough tne [loverenaimed Unddirectional and periodic, 
its oscillatory motion was not quite sinusoidal, as discussed 
in Section IV. This aberration in the flow appeared to 
manifest itself differently for high and low frequencies of 
oscillation, irrespective of the value of i. 

At low frequencies of oscillation (Fig. 8) the curve of 
BU /2 appears symmetrical about the vertical with the abso- 
lute values of acceleration and deceleration about equal. 
However, there is a long period in the cycle at which the 
valve appears to have no effect on the flow. At the hignest 
Pemerton of the cumve the valve is in the horizontal position 
with respect to the conduit axis (minimum head loss), while 
at the lowest portion of the curve the valve is vertical 
(maximum head loss). It is thus apparent that the available 
head is not significantly changed when the valve is about 
+ 30 degrees from the horizontal, causing a flat portion 
in ene Curve . 

Pa haoheaerequenciessof.oseidlations(hie, 9) oothis."null 
portion" of valve travel creates a pressure curve which is 
Unsymmetrical about the vertical characterized by a small 
acceleration and rapid deceleration. 

several attempts were made to alleviate this problem. 


First, a baffle plate was installed upstream of the valve 


50 


to decrease V. At this time vertical ridges extending 
downstream were attached around the holes in the baffle in 
order to decrease the magnitude of the deceleration. Second, 
small vertical plates were attached perpendicular to the 
valve faces so that the valve might affect the head loss in 
the near horizontal position. Finally, two plates were 
attached perpendicular to the valve faces with the same 
diameter as the valve itself to form a four-bladed butterfly 
vee. Untortunately, timeedid not permit manufacturing the 
plates to the same width or contour as the existing valve. 

i iW for this reason thatconly thosesportions of a 
Cyele whichemost nearly exhibited a sinusoidal welocity 
variation were used in evaluating the various coefficients. 
A butterfly valve system consisting of several plates is 
anticipated to alleviate this problem and yield a fairly 
samusonudal seseilwation in velocity over Ghe entire cyche. 
B. EXTENSION AND COMPARISON OF DATA WITH THOSE OBTAINED 

BY. OTRERS 

iWac deme Obbained herein are presented together with 
SMe obtained by Davenport [Ref. 36] at the University of 
Bristol and by Sarpkaya [Ref. 31] at the Aerodynamische 
Versmehsanstalt, Gottingen, msing aecirewlar eylinder of 
J imieh diameter in an oscillating wind tuamel within the 
range OI Reynolds numbers from 5000 to 9Om000. 

rarso, sdemissapparent that-thesdata of Davenportwand 
Sarpkaya constitute extensions of the present data to higher 


values of the frequency parameter, i. 


Sul 





me coma. i temomecleagrethatwctm of eemendata correlate fairly 
well with A over a large range of 4 and Reynolds numbers 
me@ardiless of whether it is obtained»with air or water in 
eeomealle orm Large tunnels 

Mined, iteis evdident stiat the eitectsmof the free—stream 
oscillations on Cy 6A) become significant only when the 
frequency of oscillation approaches the natural vortex- 
shedding frequency, i.e. when A approaches the Strouhal 
mumber of about 0.2. whlthough a theoretical explanation of 
the decrease of Cy 6A) as X} + St is beyond the scope of the 
Bmesent Study, it is not too much of a conjecture to state 
meee Che Lormmation, growth, “and motion of the vortices are 
Seeniticantly altered by the oscillations »when the shedding 
frequency is in the neighborhood of the oscillatory frequency 
of the free stream. When the period of oscillation is several 
waders Ol Moemruce larger than the periodmot vertex shedding, 
the flow behaves essentially as a juxtaposition of instan- 
Vaneous duasi-steady states and the total resistance can 
be represented accurately by the coefficients C1, (0) and 
C(A). | 

Finally, it is evident that C56) decreases sharply in 
the range of fd/V from 0.25 and 0.5. The reasons for this 
are as follows. It is a well-known fact that the drag 
acting on a cylinder undergoes oscillations with a frequency 
Wmce the fieduency of vortex shedding. In other words, as 
Ute wiiteercmes Changes 105 dircerton, up oredown, the drag 
force oscillates with an amplitude of AC. about a mean value 


D 


De 





D° Denoting the oscillation frequency of drag with 
f,, we then have f,d/V * 0.44. The same parameter for the 
ieee oscillations is f,d/V ~ 0.22 within the range of 


ei 


Reynolds numbers studied. Thus the data presented in Fig. 14 
show that Cy) is, and should be, affected in the range 
@retd7V trom 0.25 to 0.5, 1.e. as the drag—-oscillation 
frequency nearly synchronizes with the free-stream oscillation 
fmequency. Only a complete analytical investigation, if 


at all possible, can explain as to why CH OA) should decrease. 
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VI. CONCLUSIONS 


ime seoncimed onsewhwehwoeaen be drawn from this study are: 

1. The frequency parameter A» correlates fairly well the 
data obtained by various investigators regarding the resis- 
tamee coerificients. 

C. Ineephem@eangemof 0 < Ame Om05, theseresistance 
eoefficients are not materially affected by the free~stream 
oscillations. This range has not been previously investi- 
gated. All or most of the previous studies have dealt with 
eseim lations whthea zeroemean velocityged.e. with A> ©, 

per ne draeeand ineYrvia coeciitements which are functions 
of the frequency parameter A are significantly reduced in 
the range where the free stream oscillations nearly coincide 
wath the flow induced oscillations of lift amd drag. 

4, The reasons leading to the decrease in Cp 0A) and 
Cy OA) in the vicinity of A = 0.4 and A = 0.2 respectively, 
can be clarified only by a combination of additional theoret- 
ical and erinental investigations where the full Navier- 
Stokes equations are solved for an oscillatory flow about a 
circular cylinder and where the characteristics of flow in 


the near-wake of the cylinder are measured. 
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VII. RECOMMENDATIONS FOR FURTHER STUDIES 


Hhneefollowinmg are recommendatvons forsfurther=studies: 

1. Improve the butterfly valve system so that nearly 
sinusoidal oscillations may be obtained in the range of 
cae emer Cer) 

2. Undertake a study of the solution of the Navier- 
Stokes equations via finite-difference or finite-element 
methods for the range of 4» values cited above. 

3. Carry out experiments with flat plates and other 
Skiapesmotimbodies tComascertedan the: universality~of «the 


parameter i. 
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